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Results from a 26 shot ballistic range test of the Mars Science Laboratory (MSL) entry 
capsule are presented. The supersonic pitch damping properties of the MSL capsule were 
characterized between Mach 1.35 and Mach 3.5 and total angles-of-attack from 0° to 30°. 
In flight, the MSL entry capsule will utilize a radial center-of-gravity offset to produce a 
non-zero trim angle-of-attack. This offset trim angle will produce lift, enabling the capsule 
to fly a guided entry and reducing the landing footprint dimensions to within 10 km of the 
desired landing site. A lifting configuration could not be tested at the ballistic range used 
for this test as the models would swerve into the range walls, possibly damaging cameras, 
the coordinate reference system or other facility assets. Ballistic (non-lifting) data was 
extracted and will be implemented in a conservative fashion to ensure that the dynamic 
stability characteristics of the flight vehicle are bounded. A comparison between the MSL 
pitch damping results and the dynamic model of the Mars Exploration Rover capsule shows 
generally close agreement with no significant differences in damping characteristics due to 
the change in backshell geometry. Dynamic moments are also compared to the MSL reac- 
tion control system (RCS) control authority to show the controller has sufficient margin to 
easily damp any dynamic stability effects. 
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Nomenclature 


A 

Amplitude constant 

V 

Euler- Cauchy frequency 

C A 

Axial force coefficient 

t 

Amplitude growth constant 

c La 

Lift slope coefficient 

p 

Density 

Cm a 

Pitching moment slope coefficient. 


Euler angles 

Cm q 

Pitch damping coefficient Z c <$, + fr 

9(tv) d (w) 

u 

Oscillation frequency 

Cm q 

Effective (averaged) pitch damping 

Subscripts 


C nr 

Yaw damping coefficient 

CE 

Constant energy 

d 

Diameter 

eg 

Center of gravity 

I 

Moment of inertia 

o 

Peak amplitude 

M 

Mach 

q 

Pitch rate 

171 

Mass 

ref 

Reference value 

s 

Reference Area 

T 

Total 

T 

Torque 

Trim 

Trim angle of attack 

U 

Uncertainty 

xx, yy, zz 

Moment of inertia axes 

V 

Velocity- 

00 

Frees tream 

x,y,z 

Body coordinate axes 

Superscripts 


Greek 


A 

Adder 

a 

Angle of attack 

M 

Multiplier 

6 

Phase shift constant 



e 

Sine of total angle 



P 

Amplitude growth exponent 




I. Introduction 

The Mars Science Laboratory will land on Mars in 2012. The 4.5m diameter entry vehicle (EV) will 
be the largest flown to Mars, carrying the largest payload, a 900 kg rover. The MSL capsule geometry is 
similar to other recent Mars entry vehicles. Each has been based on the original 70° sphere-cone Viking 
capsules which landed in the 1970’s. Mars Pathfinder changed the backshell shape to accommodate the 
payload packaging. The Pathfinder shape was used by the Mars Exploration Rovers and Mars Phoenix, 
but MSL has again changed the backshell to accommodate the payload. Blunt bodies going back to Viking 
and earlier have found dynamic instabilities that arise as the capsules decelerate to Mach 3.5 and below. 1 
These instabilities are thought to be driven by a wake pressure lag as the capsule oscillates. 2 Therefore, 
the dynamic stability characteristics of blunt bodies can be very dependent on backshell geometry. Prior 
experience with the backshell designs of the Stardust and Genesis entry capsules 3- 4 demonstrated that small 
modifications can produce large changes to dynamic stability. As MSL has a new backshell shape, it was 
decided to obtain experimental free flight data of the MSL configuration. 

MSL will be the first vehicle since Viking to have a control system actively damping rates during entry, 
and the vehicle will be the first to maneuver in the Martian atmosphere. The EV will perform a series of 
bank maneuvers, controlling the direction of the lift vector to manage energy while flying a guided entry 
towards the desired landing site and attempting to maximize the altitude at which parachute deploy occurs. 
To properly interpret preflight simulations of the MSL entry, it is important to understand the relative 
contributions to capsule stability from its inherent dynamic stability and the corrective control authority of 
the RCS. 

Building on Stardust and Genesis experience, the Mars Exploration Rovers mission chose to determine 
pitch damping characteristics from free-flying models. Viking had used forced oscillation testing to measure 
pitch damping and Mars Pathfinder (MPF) had used that data for preflight predictions. The MPF entry 
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resulted in oscillation amplitude growth significantly exceeding preflight predictions. 5 It is believed that 
sting effects can corrupt the mechanisms (wake flow structures) that cause dynamic instabilities in blunt 
bodies and this resulted in the Viking data underpredicting the MPF dynamic instabilities. 6 A simulation 
using MER pitch damping data to reproduce the MPF entry matched much better than did the preflight 
simulation using Viking forced oscillation data. 7 Mars Science Laboratory has again returned to the ballistic 
range to ensure that the geometry changes have not introduced any new dynamic instability issues. Data 
reduction methods have improved continually as more entry vehicles seek to quantify their dynamic stability 
using ballistic range techniques. 

One area of concern for MSL is the fact that it will be trimming at a non-zero angle-of-attack. While some 
facilities can test lifting vehicles by operating at low pressures, the pressure reduction required to reduce lift 
also results in low Reynolds numbers. Free flight data was still judged to be preferable to forced oscillation 
data. MSL tested non lifting models and is applying extra conservatism to account for lifting effects. 

II. Facility and Methods 


A. Test Facility 



Figure 1. Looking uprange at the Eglin Aeroballistic Research Facility 


Freeflight dynamic stability data was collected for the MSL EV at the Aeroballistic Research Facility 
(ARF) at Eglin Air Force Base, Florida. 8 This facility is an enclosed concrete structure open to ambient 
atmosphere, near sea level conditions. The ballistic range consists of a gun room, blast chamber and in- 
strumented range. Models are propelled to test velocity in a three-inch (76mm) artillery gun. To launch 
the MSL EV models from the gun, they are held in four-petal sabots made of polycarbonate plastic. The 
gun is fired into the blast chamber where all spent propellant and blast products are collected. As the 
model/sabot package enters the blast chamber, the sabot petals peel off and are arrested by steel plates 
at the far end of the chamber. The model passes through an aperture into the instrumented range. The 
instrumented range (see Figure 1) is a 207m hallway. The walls are 3.66 x 3.66nr square for the first 69m of 
flight, then open to 4.88m x 4.88m for the remaining distance downrange. Fifty sites along the range house 
orthogonal spark-shadowgraph stations. The fields-of-view of the shadowgraph stations are 2.13m diameter 
circles. As the models fly downrange they break IR switches at each station that initiate spark light sources 
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which produce shadows of the model and shock structure on reflective screens opposite the spark sources. 
Cameras located along one hallway wall and in a pit in the floor of the range take images (shadowgraphs) 
of each screen. Beads placed on catenary wires provide surveyed reference points in each image to precisely 
locate the position and orientation of the model at each station. A chronograph is synchronized with the 
spark source discharge to mark the time each image was captured. The aerodynamic coefficients which best 
recreate the observed motion of the model are then extracted by parameter identification techniques. 9 

B. Model Geometry 

The ballistic range model tested is a scaled version of a slightly simplified MSL EV geometry. All surface 
protuberances on the flight EV were considered too small to require modeling for dynamic stability testing. 
Therefore, the MSL EV model is a body of revolution approximating the overall shape of the vehicle. The 
models tested had a overall diameter of 65.05mm and an overall length of 41.56mm. This is a 1.44% scale 
model of the flight vehicle. To get the axial eg position correct, the model was made in two pieces of different 
densities. The forebody was made of steel, while the aftbody was made of tungsten and press-fit into the 
forebody. Figure 2 shows the model geometry and assembly of the two components. Table 1 shows the 
design mass properties as given by the 3D CAD package used to design the models. The mass properties 
are based on the solid model and measured densities of the metals used to fabricate the models. Figure 3 
shows the ballistic range model and the sabot used for launching the models. 



Figure 2. Cross-section of ballistic range model (dimensions in millimeters) 


Table 1. Design Mass Properties of Eglin Ballistic Range Model 


Model Mass Properties 
(designed) 

Mass 

(g) 

I XX 

(g — mm 2 ) 

( g — mm 2 ) 

x cg 

D 


669 

2079.2 

1459.3 

0.2991 

Material Density 

Steel 

Tungsten (Billets 1-3) 

Tungsten (Billet 4) 


(. g/cm 3 ) 

( g/mm 3 ), measured 

(g/mm 3 ), measured 


0.007850 

0.01689 

0.01787 



There was no attempt to achieve a radial offset for these models. Pre-test analysis showed that lifting 
flight would result in swerving trajectories that could not be aimed to fly any significant distance down 
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the ARF without striking the hallway. Therefore, all test data is for a non-lifting geometry. Just prior to 
parachute deploy, the MSL EV jettisons a number of balance masses to move the center of gravity back to 
an axisymmetric, non-lifting configuration. From Mach 2.5 down to Mach 2.0 the MSL flight eg location is 
incrementally migrating to match the configuration tested in the ballistic range. 



Figure 3. MSL EV ballistic range model and sabot 


C. Test Matrix 

The matrix of test shots was selected to evenly populate an oscillation-amplitude/Mach space where dynamic 
instabilities have historically been observed and through which the MSL entry vehicle will be flying. The 
Mach range of interest tvas from approximately Mach 3.70 down to Mach 1.5. The angle-of-oscillation range 
varied from 0° to over 30° . The polycarbonate sabots that held the models in the gun barrel, were fabricated 
with model cavities canted at different angles relative to the axis of the gun, allowing the models to be held 
and launched at fixed initial angles-of-attack. This technique provides fairly accurate control over the initial 
oscillation amplitude. With proper sabot selection and gunpowder mass (selected from historical tables 
specific to the gun), the test matrix achieved an even data distribution over the desired range of oscillation 
amplitudes and velocities. As is typical, low speed, low angle data was difficult to obtain. Any small 
disturbances tend to grow due to dynamic instabilities and decreasing dynamic pressure as the models slow 
down. The lowest oscillation amplitude data is always obtained at the first data stations. Sabot separation 
induces some amount of oscillation, so even the first data points for 0° sabot shots often see amplitudes of 
3 — 5° or so. The completed test matrix is given in Table 2. 


An example of the shadowgraphs taken at each station are shown in Figure 4. The images were both 
taken at station 18 during shot S070325 (second of the test matrix) at Mach 2.7. The strong density gradients 
behind the bowshock cause the forebody shadow to look deformed. Therefore angle measurements are based 
on the backshell shape. The accuracy with which angles are measured in each shadowgraph is ±0.5°. 

D. Data Reduction Methods 

Two independent data reductions were performed on data collected in the ballistic range. Arrow Tech 
Associates and AerospaceComputing Incorporated are contractors that have worked closely with the Eglin 
facility for many years, developing the software for data capture as well as data reduction. Both used their 
proprietary parameter identification codes to extract static aerodynamics and pitch damping characteristics. 
AerospaceComputing authored the Comprehensive Automated Data Reduction and Analysis (CADRA 10 ) 
software used by the ARF to measure position and orientation angles from the shadowgraphs of each shot. 
Therefore, AerospaceComputing assumed the task of extracting raw position/orientation data. This data 
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Table 2. MSL EV Ballistic Range Test Matrix, March/ April 2007 


Shot No. 

Model No. 

Mass 

(M 

X C g 

D 

Iyy,zz 
(. g - cm 2 ) 

Mi 

( g — cm 2 ) 

Sabot Angle 
(deg) 

£070324 

112 

672.36 

0.299 

1474.00 

3.48 

10° 

£070325 

121 

670.14 

0.299 

1489.72 

3.67 

20° 

£070326 

124 

668.74 

0.299 

1488.31 

3.59 

O 

O 

CO 

£070327 

102 

668.64 

0.300 

1470.06 

3.62 

20° 

£070328 

109 

671.03 

0.299 

1473.50 

2.69 

0° 

£070329 

107 

663.88 

0.299 

1454.84 

2.59 

0° 

£70330 

116 

673.19 

0.299 

1477.69 

2.68 

10° 

£070331 

101 

672.32 

0.299 

1475.01 

2.67 

20° 

£070332 

126 

667.53 

0.298 

1485.83 

2.62 

CO 

o 

o 

£070333 

105 

673.46 

0.300 

1476.57 

3.52 

0° 

£070334 

110 

671.10 

0.299 

1469.94 

3.59 

0° 

£070335 

115 

671.75 

0.299 

1475.97 

3.51 

10° 

£070336 

127 

667.94 

0.298 

1490.02 

3.48 

CO 

o 

o 

£070337 

114 

669.35 

0.299 

1469.61 

2.65 

10° 

£070338 

120 

669.83 

0.299 

1471.38 

2.59 

20° 

£070339 

125 

669.21 

0.298 

1489.87 

2.46 

30° 

£070340 

108 

672.76 

0.299 

1474.48 

2.10 

0° 

£070341 

111 

671.31 

0.299 

1470.78 

2.14 

10° 

£070342 

117 

673.08 

0.299 

1476.73 

2.04 

20° 

£070343 

103 

670.88 

0.299 

1471.48 

2.10 

30° 

£070344 

106 

671.23 

0.299 

1467.64 

1.49 

0° 

£070345 

119 

671.19 

0.299 

1471.77 

1.60 

0° 

£070346 

118 

673.16 

0.300 

1477.85 

2.62 

0° 

£070347 

104 

668.64 

0.299 

1464.20 

2.63 

CO 

o 

o 

£070348 

107 (2 nd ) 

663.84 

0.299 

1454.84 

3.12 

0° 

£070349 

108 (2 nd ) 

672.77 

0.299 

1474.48 

2.85 

0° 



Hall Image 


Pit Image 


Figure 4. Shadowgraphs of shot S070325 at station 18, Approximate conditions: Mach = 2.7, a = 13°, 0 = 0° 
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was then used by both AerospaceComputing and Arrow Tech Associates for independent analyses and 
identification of aerodynamic coefficients. 

Arrow Tech Associates 11 

Arrow Tech uses the Aeroballistic Research Facility Data Analysis System (ARFDAS) 12 parameter identifi- 
cation code to extract aerodynamic coefficients from ballistic range data. For MSL and other blunt vehicles, 
ARFDAS first uses 6-degree-of-freedom (DOF) analysis of small sections of ballistic range shots, grouped 
by angle-of-attack and Mach number, to provide linear pitch damping fits over these small sections. After 
inputing the ballistic range data including time, position, attitude of the model at each station, mass prop- 
erties, and ambient atmospheric conditions, ARFDAS performs a sectional analysis to roughly map out the 
behavior of each coefficient as a function of angle-of-attack and Mach number. The full nonlinear 6-DOF 
analysis provides the final aerodynamics. From the sectional fits, an estimated table of coefficients is supplied 
to the 6-DOF program. The 6-DOF analysis uses a Maximum Likelihood Method to iteratively adjust the 
aerodynamic coefficients, finding the best fit. Adjustments required to capture nonlinear dependencies on 
angle-of-attack and Mach number are also incorporated in this analysis. ARFDAS also has the ability to 
perform the 6-DOF analysis on several similar data sets to get the best fit among them. 

AerospaceComputing Inc . 13 

The CADRA program used to scan images is also the AerospaceComputing data reduction tool used to 
extract aerodynamic coefficients. Similar to the ARFDAS analysis, AerospaceComputing uses linear theory 
to roughly map out nonlinear pitch damping curves versus Mach number and angle-of-attack. A nonlinear 
analysis follows to provide the final aerodynamics. AerospaceComputing uses a three step process. First, 
runs are divided into short segments, where velocity and oscillation amplitude may be considered constant. 
These small segments help establish aerodynamic trends with Mach number and angle-of-attack. Next, 
simultaneous linear analysis of the shots with small oscillation amplitudes helps define the nonlinear behavior 
of the models at small angles-of-attaek. Finally, nonlinear analysis is applied to all of the shots. This analysis 
involves a study to determine the functional form of the pitch damping curve, followed by the determination 
of nonlinear aerodynamic coefficient curves by fitting 6-DoF simulations through all the shots simultaneously. 
Figure 5 shows an example of a CADRA reconstruction, fit through the raw orientation data for one of the 
shots. This reconstructed trajectory uses the final extracted aerodynamic coefficients described in the next 
section. 



Figure 5. Example of CADRA2 data reduction compared to raw angle measurements for shot S070325. 
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III. Ballistic Range Results 


The 26 ballistic range shots provided a large data set from which pitch damping characteristics of the 
model were determined. At the higher velocities tested, where Viking, MPF, and MER data suggested there 
would be dynamically stable behavior, the ballistic range shots did show the capsule to be dynamically 
stable. Overall, the two data reduction efforts identified similar dynamic stability characteristics. When 
applied to trajectory simulations, both sets of dynamic aerodynamics produce behavior much like the MER 
dynamics. The actual curves, which are functions of local angle of attack, appear quite different in form. 
The curves from both data reduction efforts will be presented, followed by a discussion of these differences. 

As predicted, the ballistic range tests show the MSL capsule to be dynamically unstable at small angles 
starting between Mach 2.8 and 3.25 . This instability increases as the model decelerates (data presented 
goes down to Mach 1). The data extracted by Arrow Tech Associates is very similar to MER ballistic range 
data and found to be slightly more dynamically stable than the MER capsule. The AerospaceComputing 
results show some new features not identified in past ballistic range testing. A very rapid change in low-angle 
dynamic stability was identified as the models decelerated through a specific Mach number. Below this Mach 
number, a narrow region of dynamic instability with a large peak exists at a = 0°, with positive damping 
at all angles away from the localized peak. As the MSL EV decelerates further, the peak of the dynamically 
unstable region decreases, but extends to larger angles of attack. The AerospaceComputing analysts feel 
that this behavior may have been present in the MER data and its identification here reflects improvement in 
their data reduction analysis capabilities (augmented/guided by recent forced oscillation test data gathered 
for the Orion Crew Exploration Vehicle (CEV)). 14 In contrast, the Arrow Tech analysis was conducted in 
much the same manner as has been done for past tests, with similar results. 

It should be noted that the parameter extracted by both Arrrow Tech and AerospaceComputing is the 
combined damping coefficient, C mq + C m& . The separate damping terms due to pitch and angle-of-attack 
rates are not easy to separate in ballistic range testing and they are typically lumped together. C mq is used 
here as shorthand for the full coefficient. 

A. Static Aerodynamics 

Static aerodynamics extracted by AerospaceComputing and Arrow Tech are in very good agreement with 
earlier ballistic range test results for the MER entry capsule. 7 The MER ballistic range data is also in 
good agreement with CFD-based MER aerodynamic database 15 and old Viking data. 16 For brevity, the 
MSL static force and moment data is not presented. Refer to the earlier MER or Viking data for the static 
characteristics of 70° sphere-cone blunt bodies. 

B. Arrow Tech Associates Damping Results 

The expression for the pitch damping coefficient used in the 6-DOF analysis is: 

Cm, = + C mq 2 £ 2 + C mqi £^ + C mq m (A / Ii — M re f) + C mq m2 (Vf,; — M re f ) 2 

where 

lv? + V 2 

£ = sin(a T ) = y ~ 

and 

M ref = 2.5 (3) 

M re f is an arbitrary constant used in the data reduction. A table of G mqinitial coefficients are determined 
from linear fits using sections of individual shots. This table is an initial estimate of the nonlinear pitch 
damping coefficients. The table of initial coefficients is then corrected to better resolve Mach number and 
angle-of-attack variations by fitting 6-DoF trajectories through the data points of multiple, complete ballistic 
range shots. The final table of coefficients, listed in Table 3, includes the extra terms in Equation 1 and 
reflects the identified dependencies on angle of attack and Mach number. 


( 1 ) 

( 2 ) 
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Table 3. Arrow Tech Associates supersonic Cm, 


results, Xcg/d, = 0.30 


Mach 


a 

0 

1.35 

1.5 

1.75 

2.0 

2.25 

2.5 

3.0 

3.5 

4.0 

0.9878 

0.9093 

0 . 7246 

0.5652 

0.4442 

0.3231 

0.1098 

-0.0900 

-0.2898 

2 

0.8882 

0.8127 

0 . 6402 

0.4914 

0.3753 

0.2592 

0.0441 

-0.1422 

-0.3285 

4 

0.6438 

0.5765 

0.4321 

0.3078 

0.2052 

0 . 1026 

-0.1107 

-0.2655 

-0.4212 

6 

0.3650 

0.3066 

0 . 1950 

0.0990 

0.0113 

-0.0765 

-0.2700 

-0.3933 

-0.5175 

8 

0.1457 

0.0945 

0 . 0087 

-0.0648 

-0.1409 

-0.2169 

-0.3789 

-0.4806 

-0.5823 

10 

0.0153 

-0.0315 

-0.1023 

-0.1629 

-0.2318 

-0.3006 

-0.4329 

-0.5238 

-0.6147 

12 

-0.0456 

-0.0903 

-0 . 1542 

-0.2088 

-0.2745 

-0.3402 

-0.4536 

-0.5400 

-0.6273 

14 

-0.0693 

-0.1134 

-0.1741 

-0.2259 

-0.2903 

-0.3546 

-0.4590 

-0.5445 

-0.6309 

16 

-0.0757 

-0.1197 

-0.1799 

-0.2313 

-0.2952 

-0.3591 

-0.4599 

-0.5454 

-0.6309 

18 

-0.0771 

-0.1208 

-0.1809 

-0.2322 

-0.2961 

-0.3600 

-0.4608 

-0.5463 

-0.6318 

20 

-0.0779 

-0.1218 

-0.1814 

-0.2322 

-0.2966 

-0.3609 

-0.4608 

-0.5463 

-0.6318 

25 

-0.0779 

-0.1218 

-0.1814 

-0.2322 

-0.2966 

-0.3609 

-0.4608 

-0.5463 

-0.6318 

30 

-0.0779 

-0.1218 

-0.1814 

-0.2322 

-0.2966 

-0.3609 

-0.4608 

-0.5463 

-0.6318 


The pitch damping coefficient, C mq , is plotted in Figure 6 for a range of Mach numbers from Mach 
3.5 down to Mach 1.4. These curves, interpolated from the final coefficient matrix in Table 3, show the 
capsule to be dynamically unstable (C mq > 0) at low angles-of-attack. The magnitude of instibility is seen 
to increase as the MER capsule decelerates towards Mach 1.4. At higher angles-of-attack, the capsule is seen 
to be stable but becomes less so with decreasing Mach number. The Mach variation is roughly a linear shift 
of the entire pitch damping curve. Figure 7 shows plots of total angle-of-attack and Mach number versus 
downrange position for three representative shots. The initial angles-of-attack were similar for each shot, 
so the difference in amplitude growth between them is due to the dynamic stability variation with Mach 
number. The slower the initial velocity, the greater the amplitude growth downrange. The pitch damping 
curves and behaviors observed are generally similar to the Arrow Tech Associates MER data. 



a (deg) 


Figure 6. Arrow Tech Associates supersonic pitch damping coefficient curves, C mq (xcg/d = 0.30 ,z cg /d = 0.0) 
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o 

CO 


Downrange (m) 


Figure 7. Example of amplitude growth with initial Mach number 


C. AerospaceComputing Results 

The AerospaceComputing pitch damping coefficient data is plotted in Figure 8 for a range of Mach numbers 
versus angle-of-attack. The functional form describing these curves is: 


A/qo 2.8, — TTldX 


(C mqA + C rriq2 M 00 + C mq3 Ml + C mqA Ml sin 2 a) , (C mqA + 


Mcx> > 2.8, Cm q — Cm, i5 + 


c. 


m q , 6 

Ml, 


( 4 ) 

( 5 ) 


Coefficients for these functions were determined by multiple fits of the ballistic range shots with 6-DOF 
nonlinear analysis. These coefficients are 


C mqA = 1.5733 
C mqt2 = -0.9797 
C mqi3 = 0.4985 
C mqA = -30.9 
C mq , 5 = -0.47 
C mq , 6 = 0.16 

The curves show interesting behavior with angle of attack and Mach number. The AerospaceComputing 
results indicate that there is a distinct Mach number as the MSL EV decelerates, M =2.8, where a dynamic 
instability arises. The peak of this instability at a = 0° is very large, but only extends about two degrees 
from the peak. At lower Mach numbers, this peak decreases, but extends to larger angles-of-attaek. This 
instantaneous onset of dynamic stability with a very narrow, but steep peak has not been identified in other 
ballistic range data analysis (MER, Genesis, Stardust, etc.). However, features like this have recently been 
identified in forced oscillation testing of the CEV (Orion) command module. This is thought to be reflective 
of the advancements in data reduction technique and the increasing knowledge base of the data analysts, 
rather than a distinct behavior exclusive to the MSL entry vehicle. While the functional form differs greatly 
from the Arrow Tech Associates results, the integrated effects of both sets of data, when used in simulation, 
replicate the measured data sets very well. A comparison of the integrated effects is presented next. 
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Figure 8. AerospaceComputing supersonic pitch damping coefficient curves, C mq (x cg /d = 0.30 ,z cg /d = 0.0) 


D. Discussion of Arrow Tech and AerospaceComputing Results 

The functional forms of AerospaceComputing and Arrow Tech Associates results are quite different. While 
the Arrow Tech results are very similar to MER results, the form of the AerospaceComputing damping curves 
are very different than the measured dynamic stability curves of past Mars and Earth return missions. To 
compare them, the mean value theorem is used to determine the integrated effect of dynamic damping for 
a given oscillation amplitude. This integrated effect allows the user to replace a very nonlinear curve with 
a constant value. At a given oscillation amplitude, the constant value will produce oscillation growth or 
decay identical to the nonlinear curve. This provides a useful way to compare two very different nonlinear 
damping curves and assess how they will drive oscillatory behavior in flight. The application of the mean 
value theorem to express damping as function of oscillation amplitude was well documented by Redd et 
al. 17 The exact application used here is described by Schoenenberger and Queen. 18 The effective dynamic 
damping for a decelerating blunt body in a constant density flowfield is given by 


Cm „ = 


Jo“ (P n 


C„ 


Jo 


y 2 dt 


(7) 


Figure 9 shows the effective damping versus oscillation amplitude for the AerospaceComputing and Arrow 
Tech results. At Mach=2.8 and lower, the AerospaceComputing results show the dynamic stability to be 
much more unstable at small angles of attack. However, the AerospaceComputing curves collapse to behave 
much like the Arrow Tech data when oscillating at amplitudes of 10° or greater, and the effective damping is 
almost independent of Mach number. The low angle-of-attaek damping trends toward becoming dynamically 
stable as the capsule decelerates, following the discrete jump at Mach=2.8. When used in simulation, the 
AerospaceComputing and Arrow Tech curves will add energy over approximately the same small angle range 
and remove energy over the same large angle range. The large positive C mq values seen at the small angles 
of the AerospaceComputing curves will cause the capsule to grow to neutral energy oscillation amplitudes 
much faster than the Arrow Tech curves. This is the chief difference between results from these two data 
reduction efforts. The results are in closest agreement at moderate oscillation amplitudes ( a = 10° to 20°) 
near Mach 2.5. The majority of ballistic range shots passed through this oscillation-amplitude and Mach 
space and the curves are therefore anchored by the most number of data points in this space. 

The oscillation amplitudes where C rHq crosses zero is the point where the dynamic damping does not add 
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Figure 9. Effective damping for several Mach numbers, (x cg /d = 0.30. z cg /d = 0.0) 


to or subtract from the oscillatory energy of a body in free flight. A plot of these constant energy amplitudes 
for both sets of data is shown in Figure 10. Note the crossing point where the two sets of data agree is 
at Mach 2.4 and an oscillation amplitude of about 15°. The discrete jump of the AerospaceComputing 
data at Mach=2.8 is also obvious. Above that point, the AerospaceComputing data indicate that capsule 
oscillations are losing energy at all amplitudes. The Arrow Tech data indicate a similar falloff in constant- 
energy amplitude with increasing Mach number. The capsule must be at speeds above Mach 3.27 before the 
dynamic term is stable for all angles-of-attack. At low speeds, the Arrow Tech data indicates the capsule 
must be oscillating about amplitudes much greater than those of the AerospaceComputing data to achieve a 
constant energy oscillation. However, as shown in Figure 8, the AerospaceComputing damping model causes 
much faster amplitude growth at low angles, up to the constant energy amplitude. In the ballistic range, 
the Mach number is changing very rapidly and the flight time is so short that oscillation amplitude growth 
or decay generally never reaches equilibrium during a shot. The determination of the exact functional form 
with this test technique may not be possible. In general, where there is the most data to anchor the curves, 
these two analysis efforts are in closest agreement. For use in flight simulation it is conservative to use 
the curves which produce the most destabilizing behavior. Since the MSL vehicle is controlled, it will be 
oscillating about small amplitudes. Therefore, the AerospaceComputing data which causes rapid growth up 
to moderate amplitudes is the more conservative set of data. Above Mach 2.8, it is more conservative to use 
the Arrow Tech data which still has regions of dynamic instability. 

The MER aerodatabase constant energy curve was also calculated and included in Figure 10. As was 
shown in the Figure 6, the MER and MSL Arrow Tech data is very similar in form, differing roughly by a 
constant value. The MER data reaches a constant energy point at larger angles of attack and becomes fully 
stable at Mach 3.5. To be conservative, the MER data will be used in the MSL aerodynamic database at 
Mach numbers greater than 2.8. 

The constant energy amplitude point is where dynamic stability adds no energy to the system. However, 
as the MSL capsule is decelerating, oscillation amplitude still increases as the drop in dynamic pressure re- 
duces the restoring static moments acting on the vehicle. The oscillation amplitude growth for a decelerating 
body in a constant density flowfield can modeled by the Euler-Cauchy equation as described in two papers 
by Schoenenberger and Queen . 18,19 The solution of the E-C equation is of the form: 

a = cosfVlnf + S) (8) 
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Figure 10. Constant energy amplitude versus Mach number. ( x cg /d = 0.30, z cg /d = 0.0) 


Where S is phase shift angle, v is a frequency term, and fi is a damping term. The oscillation amplitude 
of a decelerating capsule grows proportional to time raised to the damping term, /i. This term is: 


1 C La j md 2 ^ — 
2C a + AIC a m " 


( 9 ) 


While, C mq = 0 results in constant oscillatory energy, it does not result in constant amplitude oscillations. 
To achieve constant oscillation amplitude, the dynamic damping must be sufficient to negate the amplitude 
growth due to deceleration and lift generated as the capsule oscillates. Setting /i equal to zero the C mq value 
required for constant oscillations is found to be 


C' rn i ,a=constant ^ ^ TYld 2 (^ 9 ) 

The constant amplitude damping is equal to about —0.27 for the MSL entry vehicle. 

When the effective damping satisfies n= 0, the capsule will see no oscillation amplitude growth. Equation 
9 shows that the effective damping must be negative to counteract the constant 1/2 term which produces 
natural oscillation growth due to decreasing dynamic pressure as the EV decelerates, and the Cl 0 term 
which causes growth from a coupling between heaving due to lift and the oscillations of the capsule. As the 
body oscillates, the lift force (itself dominated by the pointing of the axial force vector) produces heaving 
motion that adds an induced component to the angle-of-attack rate. Through this induced rate from the 
heaving motion, the static aerodynamics alone can affect amplitude growth. 


E. Comparison to the Mars Exploration Rover Entry Capsule 

It is useful to compare the MSL ballistic range data to the MER ballistic range data used in the MER 
flight aerodatabase. As MER and MSL have essentially identical forebody shapes, this comparison shows 
the impact of changing backshell geometry on the dynamic stability of the vehicle. Figure 11 shows the 
effective damping of the two capsules as identified by Arrow Tech Associates. Arrow Tech’s MER ballistic 
range results were used in the MER aerodynamic database. The damping functional forms are essentially 
identical, but the MSL capsule is slightly more stable than MER. The two sets of curves differ by a constant 
offset roughly equal to a change in Mach number of about 0.2 to 0.3. 


13 of 19 


American Institute of Aeronautics and Astronautics 


The functional forms were not specifically selected to be of identical form. The similarity of the curves 
is likely a result of commonality in methods and personnel used by Arrow Tech Associates to extract the 
dynamic damping coefficients. It is not likely than any physical phenomena can be inferred from the similar 
results, but this consistency in data reduction analysis adds confidence that the overall MSL dynamic stability 
characteristics are similar to MER and slightly more stable. 



Figure 11. Comparison of MER and MSL effective damping for several Mach numbers, ( x cg /d = 0.30, z cg /d = 0.0) 

The functional form identified by AerospaceComputing for the MER ballistic range data was similar 
to the MER Arrow Tech results. The new functional form identified for the MSL shape is believed to 
reflect improvements in dynamic stability modeling overall (guided by recent improvements in blunt-body 
forced oscillation tests) and not a distinct change in the dynamic stability characteristics of the vehicle. A 
comparison using MER and MSL pitch damping models to replicate the MSL data points was performed to 
assess their differences. These comparisons showed the integrated effect of the pitch damping curves. While 
local discrepancies, especially at low oscillation amplitudes (where there are the fewest data points anchoring 
the models) were observed, in general both models fit the MSL data well. The comparison suggests only small 
differences in the overall dynamic stability characteristics of the MER and MSL vehicles. The differences in 
the pitch damping functional forms were not driven by any significant changes to dynamic stability because 
of backshell geometry changes. 


F. Uncertainties 

The uncertainties used for the MSL database are based on those used for MER, but their implementation 
has been modified to produce more conservative dispersions across the Mach and angle-of-attack space. The 
pitch and yaw damping dispersion model in the supersonic regime is as follows: 


C„ 


q,di , 


= (i + 0(c, 


+ 0.5) - 0.5 + U t 


( 11 ) 


n — n 

Hr, disperse TTLq , disperse 

Where the multiplier and adder dispersion values are: 

t/ g M = [-50% to + 150%] = [-0.5, +1.5] 


( 12 ) 


( 13 ) 
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U* = [0.00 to ± 0.10] (14) 

The 0.5 terms in the dispersion model are used to keep negative regions of the pitch damping curves 
from becoming more negative when a positive dispersion is applied. The curves are shifted up so the pitch 
damping is zero or positive, the multiplier is applied, and the curves are shifted back down by the same 
increment. The adder is then applied to disperse the asymptote. The adder value of 0.10 was used in the 
MER database as a dispersion on the Mach variation of the pitch damping curve. A value of 0.1 corresponded 
to a change in Mach number of approximately 0.5. The multiplier was selected to disperse the peak dynamic 
instability. As shown in the two MSL data reduction efforts, the pitch damping at small angles-of-attack is 
the region with the greatest uncertainty. The multiplier was chosen to be rather large to account for this 
uncertainty. To be conservative, the multiplier is also biased in the destabilizing direction. The multiplier is 
dispersed using a Gaussian distribution, while the adder is distributed uniformly. 

Another issue regarding uncertainties is that the MSL capsule is trimming at a non-zero angle-of-attack 
over much of the Mach range covered by this test. At Mach 2.5 the capsule is trimming near —20°, at 
which time six entry balance masses (~ 29 kg each) are jettisoned in sequence, smoothly changing the trim 
angle to 0° just before parachute deploy at Mach 2.0. If the unaltered pitch damping model were used 
in simulation, the capsule would experience only positive damping at the large trim angles as the capsule 
would not see the peak instabilities at a = 0°. Testing with the correct eg offset would be required to 
quantify how the damping characteristics change with trim angle. For MSL simulations, the non-lifting data 
is implemented in the most conservative manner possible; the peak dynamic instability is shifted to the 
current trim angle-of-attack (accomplished by shifting the angle-of-attack variable used to call the damping 
model) . This modification conservatively accounts for possible radial eg offset effects and the accompanying 
non-zero trim angle. 

G. Comparison to Reaction Control System (RCS) Authority 

To better understand the magnitude of the MSL dynamic instabilities in the context of controlled flight, the 
dimensional torques caused by the dynamic damping coefficients are compared to the control authority of the 
reaction control system. The MSL controller uses four pairs of reaction control jets to perform maneuvers 
and damp any unwanted rates. Figure 12 shows the reaction control system on the MSL backshell. The jets 
are scarfed with the nozzle exits approximately flush with the backshell. The jet pairs are canted to provide 
components of roll, pitch and yaw. The controller selects the proper combination to generate torques about 
a desired axis. Each jet generates approximately 290 N (65 lb). The reaction control system is capable of 
generating roll and yaw torques of ±900 (jets 3,4,7 and 8, or 1,2,5 and 6) and ±600 Nm (jets 5,6,7 and 8) 
respectively. The RCS pitch authority is ±1200 (jets 3,4,5 and 6) and —1000 Nm (jets 1,2,7 and 8). The 
asymmetry is due to the radial eg offset in the pitch plane. 



Figure 12. MSL RCS jet configuration, jets are labeled with MSL numbering convention 
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To assess the dynamic torques, the damping coefficients are combined with the rates associated with a 
range of oscillation amplitudes for several points along a nominal reference trajectory. At a given Mach 
number and dynamic pressure, the static stability can be used to calculate the peak rotational rate the EV 
will see as it oscillates at the specified amplitude. A sinusoidal oscillation is assumed. The effective peak 
constant torques and the peak instantaneous torques are plotted here against the nominal MSL RCS pitch 
torques. 

The peak torque is the product of the peak pitch damping value ( C mq at arnm ) and the peak rate (over 
an oscillation) which also occurs at the trim angle-of-attack. This the largest instantaneous dynamic torque 
that the capsule will see for a given oscillation amplitude. 

Tpeak = (ct — 0 ) • Ap ea fc(o 0 ) (lb) 

It is useful to compare the peak destabilizing torques to the RCS torques, but note that away from the peak 
rate are regions where the pitch damping curve is less positive and at large angles, negative. When the EV 
is at angles away from the angle of peak pitch damping, the controller will be fighting much smaller torques, 
or helped by regions of positive damping. As shown in Equation 7, the damping can also be expressed as a 
constant value for a given oscillation amplitude. The effective damping is a constant value that can replace 
a nonlinear pitch damping curve for a given oscillation amplitude. The effective peak torque is simply 
the effective damping for a given oscillation amplitude, multiplied by the peak rate for that amplitude. 
Comparing effective torques against the controller authority shows the capacity of the controller to negate 
any undamping produced by dynamic instabilities. 


T effective — ^rn q (o G ) * OLpeak(&o) 


(16) 


The peak rate for a given oscillation amplitude can be estimated by modeling the EV pitching oscillations 
as a simple harmonic oscillator: 


.. pV^Sd 2 (n , ^ ^ . pV^Sd^ _ n 
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For constant coefficients, the solution to Equation 17 is: 


(17) 
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where a 0 is the oscillation amplitude and the frequency of oscillation is 
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Compared to static moments, the damping moments are small for a single oscillation cycle and the 
amplitude can be approximated as constant. The peak rate is then 


apeak = a 0 wcos((wt + d) 


0 ) = a o L0 = OtopiVooi-CmJ = (^-)' 


21, 


yy 


(20) 


The peak rate is directly proportional to the oscillation amplitude, the square root of the dynamic 
pressure and the square root of the static stability. Here the dynamic pressure is broken out to highlight the 
dependency of the peak rate due to density and velocity as well. 

To demonstrate the dimensional torques that dynamic stability terms produce in flight, a typical reference 
trajectory and the allocated mass properties are used to predict the peak rates for a range of oscillation 
amplitudes. The local peak torques and effective torques integrated over a full oscillation are then presented 
and compared to the RCS pitch authority of the MSL EV. The allocated mass and pitching moment of 
inertia of the MSL EV are 


m = 3067.0 kg 


I yy = 3352.62 kg — m 2 
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Figure 13 shows the peak torques that the MSL EV will see at key points along a reference trajectory. 
These torques are calculated using the a = 0° C mq data in the MSL database (AerospaceComputing ballistic 
range data for the Mach 2.0 and 2.8 curves) multiplied by the peak rate that the EV will see for a given 
oscillation amplitude, given by Equation 19. Note that the hypersonic pitch damping value is negative, so any 
rates are naturally damped at the point of peak dynamic pressure and any point above Mach 5.0. Mach 2.8 
is the point of peak dynamic instability in the aerodatabase and therefore produces the greatest peak adverse 
torque for a given rate. The +3er dispersed curve does exceed the RCS control authority, but the EV must 
be oscillating ±25° or greater about the trim point for this to occur. Recall that the pitch damping curve 
is very nonlinear with angle-of-attack at the supersonic Mach numbers, and even this large peak adverse 
torque will be mitigated significantly as the EV swings through angles-of-attack with positive damping. At 
Mach 2.0, the peak adverse torque is slightly smaller than at Mach 2.8, and the dynamic pressure has also 
decreased. This results in a significantly reduced peak torque at parachute deploy conditions. Note each 
curve in Figure 13 is calculated up to an oscillation amplitude of 30° . The a 0 = 30° peak rates at each flight 
condition (the rate at the end of each line in Figure 13b) give an idea of how the combined effects of the 
dynamic pressure and static stability vary along the trajectory, changing the influence of the pitch damping 
data. 


Peak Dynamic Pressure (Mach 18) 

Mach 2.8 




Peak Rate (deg/s) 

(b) Peak damping vs. peak rate 


Figure 13. Peak damping torques (at a = ctTrim) based on AerospaceComputing BR data at several points 
along reference trajectory 


While it is important to note the peak possible torques that the dynamic stability characteristics can 
impart on the vehicle, the integrated effect is a more representative indicator of the undamping effects a 
controller will have to overcome in practice. The mean value theorem (Equation 7) shows that the nonlinear 
pitch damping curve may be replaced by a constant value for a given oscillation amplitude. Likewise, the 
product of the damping coefficient and instantaneous angular rate can be integrated over time to express 
the dynamic torques as a constant value for a given oscillation amplitude. The effective damping torques 
can be thought of as a direct increase or reduction of the RCS control torques, directly proportional to 
the instantaneous rate. These damping effects vary with oscillation amplitude (which drives the peak rate 
and effective damping) and Mach number (which also changes the aerodynamic damping). Figure 14 shows 
effective torques versus oscillation amplitude and the peak rates that occur for those amplitudes. Again, in 
Figure 14b the curves were generated for oscillation amplitudes up to 30° . The variation in peak torques for 
the different curves is due to changes in dynamic stability and the variation in peak rates for the different 
curves is due to changes in static stability and dynamic pressure along the entry trajectory. 

Note that the +3er peak adverse torque at Mach 2.8 reaches a peak of approximately 75 Nm at a 0 = 4°. 
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For small angle oscillations, the region of dynamic instability dominates the effective damping. The effective 
damping torques grow nearly linearly due to the proportional increase in rates as amplitude increases. 
For larger oscillation amplitudes, the negative C mq values at larger angles of attack provide more positive 
damping, resulting in a gradually decreasing effective damping for any reasonable oscillation the MSL EV 
might see. At hypersonic speeds (represented here by the peak dynamic pressure curves), the damping is 
dominated by the forebody only and is conservatively modeled in the aerodatabase as a constant value. 
Therefore, the effective and peak torques are identical. 
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(a) Effective damping vs. amplitude 
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Figure 14. Effective damping torques based on AerospaceComputing BR data at several points along reference 
trajectory 


IV. Conclusion 

A ballistic range test was performed to determine the dynamic stability model for the Mars Science 
Laboratory entry vehicle. Two independent data reduction efforts identified nonlinear models for pitch 
damping, C mq , as functions of angle-of-attack (a from 0° to 30°) and Mach number (Mach 1.5 to 3.5). 
While the functional forms identified by the two efforts look very different, they predict similar levels of 
dynamic stability when used in simulation. The agreement is closest in the Mach and angle-of-attack space 
anchored by the largest number of ballistic range data points. The largest disagreement occurs at small 
angles. This is not unexpected as it is very difficult to get small oscillation data in the ballistic range. At 
larger oscillation amplitudes, the effective damping of the two models are in very close agreement. The MSL 
pitch damping was found to be very similar to the Mars Exploration Rover pitch damping characteristics 
(determined with a similar ballistic range test). MSL may be slightly more dynamically stable than MER 
and no significant effects due to the change in backshell geometry were identified. 

Lifting models cannot be safely tested at the Eglin AFB ballistic range, so all models tested were bodies 
of revolution with no eg offset. To account for lifting effects when the data is used in simulation, the peak 
instabilities are shifted to the current trim angle-of-attack. This approach is used for the MSL aerodatabase 
to be conservative. Further testing would be required to assess the true variation of pitch damping with 
radial eg offset. 

The destabilizing effects of the pitch damping data were compared to the reaction control system of the 
MSL EV. Peak dynamic torques should never approach the torques generated by the RCS system along any 
part of the trajectory and the effective damping should act only as a slight perturbance to be easily countered 
by the controller. As the controller is continuously damping any excessive rates, the dynamic instabilities 
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can not act on the vehicle uncorrected. No large oscillation amplitudes or rates are expected. 
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